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Abstract. In this paper we propose a particular description of meson 
spectroscopy, with emphasis in heavy bound states like charmonia and 
bottomonia, after working on the main aspects of the construction of an effective 
potential model. We use the prerogatives from "soft QCD" to determine the 
effective potential terms, establishing the asymptotic Coulomb term from one 
gluon exchange approximation. At the same time, a linear confinement term 
is introduced in agreement with QCD and phenomenological prescription. The 
main aspect of this work is the simplification in the calculation, consequence of a 
precise and simplified description of the kinetic term of the Ifamiltonian. With this 
proposition we perform the calculations of mass spectroscopy for charmonium and 
bottomonium mesons and we discuss the real physical possibilities of developing 
a generalized potential model, its possible advantages relative to experimental 
parameterization and complexity in numerical calculations. 
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1. Introduction 

Since the discovery of QCD (Quantum Chromodynamics), there have been remarkable 
technical achievements in perturbative calculations applied to hadrons. However, it 
is difficult to use QCD directly to compute hadronic properties. In this context, 
phenomenological potential models have provided extremely satisfactory results 
in describing ordinary hadrons, more specifically quark - antiquark bound states 
(mesons). 

The pioneer work proposing a potential model inspired in QCD was the De Riijula, 
Georgi and Glashow model [T]. Since then, the phenomenological prerogatives to 
characterize potential models were based on: (i) hadrons are in a color singlet; (ii) at 
short distances, QCD interactions obey asymptotic freedom and one gluon exchange 
is dominant with a limited coupling constant (potential will be of Coulomb type) ; {Hi) 
at long distances, the quarks are bound by a flavor independent long range potential, 
obtained phenomenologically; (iv) quarks in a bound state move relativistically and 
the approximation to first order in u^\c^ is normally used. 

Among several effective models worked out until now [H[2l[3ll4l[5|[6l.[7l[8l.[9lfT0l[TT]. 
the less complex alternative involves the Schrodinger formalism, that can also be 
derived from the Salpeter description of a bound state ([T]): 

[E-H,- hS{p, P)^ j d^p'i^tl'.Uip - p')4>{p', P)72°AJ 

-Ar7?t/(P-P')0(p',^)72%"}, (1) 

where 




(2) 

4>{P,P)= J dp''^{p\p,P), (3) 

and Hi is the Dirac Hamiltonian. 

From Hi we obtain the kinetic approximation term T**" and with the following 
approximations @ the perturbative form of Schrodinger equation (O is obtained: 

r A+[/(p - p')~4>{p')Kt U{p - p'Mp') 

[A~U{p~p')^{p')Ar ^ 



■ 2 

^ ^ rpkin 



E 



ip') = - / d'p' Uip - p')4>ip')- (5) 



The physical prerogatives impose an effective Hamiltonian only with vector and 
scalar Lorentz components. It includes all Breit - Fermi interaction terms [4l [T2l [T3] . 
consequently, all spin effect terms. 

2. The Hamiltonian kinetic term 

Since the pioneer papers of Stanley and Robson [14] and Godfrey and Isgur [6] the 
kinetic term effects of a quark - antiquark bound state are observed with emphasis 
to relativistic effects, although treating heavy quark systems. Keeping in mind the 
perturbative treatment of relativistic corrections, despite its effects making the wave 
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functions result more complicated [51 1151 1161 117j , we determine a kinetic energy term 
approximation that describes the mean behavior for low energy states. We adopt a 
nonrelativistic approximation based on Martin inequality ^ as suggested in [T71 [15] : 



M p 



2 



Jp2 + < — + -i— + . (6) 



For a given M, the equality in ^ is obtained at the point po = \/ M'^ — m? in 
the momentum space. The curves defined by the two sides of equation ([6|) are tangent 
at Po, suggesting = +Pq^. 

Since Martin's condition is an operator relation, that inequality will continue valid 
for mean values of a state and for the average of mean values over a set of states with 
arbitrary M and po. In \l7j the choice of po^ = < p^ > was suggested by Jaczko and 
Duran for a set of states, with a reasonably exact nonrelativistic approximation. The 
equation ^ can be rewritten in the following way: 



with 



M = v/< P^ > (8) 

The physical content of this result can be illustrated through the following direct 
expansion around Pq^ 



y/j^ + m^ = yjp^ - po^ + AP 

= Af + ^l^-(^— ^ + .... (9) 

2M UP ^ ' 

We observe that an approximately correct choice of the parameter will provide 
an excellent result for mass spectroscopy of heavy quarkonium systems. In [17) a linear 
approximation in p^ with a translational term e was adopted and the parameter po^ 
is suggested and calculated to be the average of the mean square values of momenta 
over the first four energy states. 

The kinetic energy of each particle in the bound state: 



T'="(e.act) = VV+^ (10) 

can be approximated by an expression which is linear in p^ as: 



,2 „ 2 



T^'\,,.) - V^^FTm^ + ^" ^ . (11) 

2VPo + TO 

Although this parameterization has the same structure as that adopted by Jaczko 
and Duran in [T7|, we followed a different procedure to calculate the values of po'^. 
The values of po^ are determined as to minimize the mean square deviation (MSD) 
of the linear approximation PT|) from the exact relativistic kinetic term PH)) . We can 
visualize the results of the MSD procedure in figures l.(a) and l.(b), in the interval 
GeV^ < p^ < 10 GeV^, depending on po^ for different values of mass interesting in 
the study of heavy quark systems (from GeV to 5 GeV) . 

The advantage of such procedure is that the obtained values of po^ are 
independent from any energy state used as input data. Consequently, they can be 
applied to calculating the whole mass spectrum of a bound state while in [17] the 
calculation is not extended beyond the first four energy states of a bound state. 
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Figure 1. (a) Mean square deviation of exact and approximated kinetic term 
for masses between and 2.0 GeV. 
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Figure 1. (b) Mean square deviation of exact and approximated kinetic term 
for masses between 2.5 and 5 GeV. 



A simple kinetic approximation for heavy meson mass spectroscopy 



5 



3.0 




1.0 I ■ 1 ' 1 ' 1 ' 1 ' 1 ' 1 

1 2 3 4 5 6 

p2 (GeV2) 



Figure 2. (a) Linear approximation and exact term versus with Pq = 
2.8 GeV^, for m = 1.2 GeV. 
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Figure 2. (b) Linear approximation and exact term versus witli pg = 
2.9 GeV^, for m = 4.2 GeV. 
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Expression pip will be extremely accurate in the region of kinetic energy where 
the probability of finding a quark in a bound state is significant. This region 
corresponds to an intrinsic energy of each particle < 10 GeV^ in the center of 
momentum frame. 

We note that the value of the parameter po^ is approximately 2.9 GeV^ for a 
particle with mass between 4.0 GeV and 5.0 GeV, as is the case of the bottom quark. 
In the case of the charm quark, with mass between 1.0 GeV and 2.0 GeV, the value of 
the parameter po^ is approximately 2.8 GeV^. 

Figures 2. (a) and 2.(b) show the exact kinetic term and its linear approximation 
for these two cases, using typical values of masses in the regions mentioned above and 
for the main region of intrinsic momenta. 

The procedure adopted here avoids sistematic errors induced by simultaneous 
evoluation of the parameter pq^ and the effective quark masses. 

Other treatments of the kinetic nonlocal terms have been already proposed 
[51 [TTJ [ini [T71 [12] , however, none of them provides so simple subsequent calculations. 



3. The interacting potential 

Our propositions will follow some prescriptions that were already developed in [6l [^ . 
We start with the Coulomb plus Linear form of the interacting potential without 
perturbative corrections. 

We know that QCD provides expression for the Coulomb coupling constant 
through a first order approximation following Feynman diagrams [20] : 

a,{k^) = (12) 

(33 -2n/) In (II) 

The scale parameter A lies in the interval between 0.15 GeV and 0.30 GeV when 
the flavor number (jif) is 4 or 5 20J. Observing the behavior of Q;s(fc^), it approximates 
to a constant value in the region of k where 5 GeV < fc < 20 GeV. However, when 
k — > A, the first order perturbative treatment is not possible since (jl2p diverges 
indicating the confinement. 

Godfrey and Isgur [6] proposed a parameterization of such behavior for 5 GeV < 
fc < 20 GeV in the following way: 

w,^ r 

as(fc^) = V'amexp — =4> (r) = ^ a„i / exp (-.x^)(ia;. (13) 

A simpler approximation can be suggested in the treatment of heavy quark 
systems: 

as(r) « as [1 - exp (-wr)] . (14) 

This parameterization correctly describes the asymptotic behavior of the coupling 
constant. Consequently we may assume it is correct just for the asymptotic region 
where the perturbative treatment may be applicable. 

The confining term will turn up in agreement with experimental data [7] that 
indicate a linear behavior. In the interaction region 5 GeV < fc < 20 GeV just the first 
term of the summation (|13p is sufficient to describe the interaction potential (|15p , and 
in function of the parameter w we can reproduce this behavior through the equation 
(|14p . This procedure is illustrated in figure [3| 

4 a 

Uir)^--—[l-exp{~wr)] + ar + b. (15) 
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Figure 3. Coulomb coupling constant parameterization in the use of 
approximation II14II . 



Obviously, expression (fT5|) just describes the pure terms of the potential (C/(r) = 
Scalar(r) + Vector(r)), but in our calculation all interacting elements with orbital 
angular momentum, spin and nonlocal term dependencies were comprised following 
the Gromes prescription [H [71 [21] for the study of the Breit - Fermi Hamiltonian \12\ ■ 

It is also important to remember that a quark bound state is described as an 
effective degree of freedom "smeared out" by gluons and quark - antiquark pairs in 
such a way that some authors propose a density function of quarks in configuration 
space describing such potential delocalization effect [SI [51 [TT]. The delocalization 
effect, in a simplified description like this one, is already achieved with the parameter 
w of the equation (|14p . 

A finite width effect parameterization could be included. It generates a 
mathematical mixing with the delocalization effect. For simplicity, we develop the 
calculations just with one parameter (w) describing both effects mixed. 

4. Methods and calculations 

Considering that our prescription involves the study of a perturbative Hamiltonian 
that comprises all relativistic corrections appearing in Gromes prescription [H [7], 
we use simple calculations in opposition to more extensive methods like orthogonal 
collocation ,15j or some not so exact variational method. We find the parameters 
choosing a simplification adopting the mean square minimization of experimental 
masses of eight resonances (four charmonia and four bottomonia) that are excellent 
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input data 

Adopting the wave function ij{r) = Rni(r)Yi„i{f), the radial equation for it is 
written like 

(r'=™ + Us{r) + U,{r) - Eni)Rni{r) = 0. (16) 

The terms Us{r) and Uv{r) involve all relativistic corrections following the Gromes 
prescription [31[7]. It permits us to obtain the eigenvalues Eni for all quantum numbers 
of a quark - antiquark system. 

The proper choice of the radial functions depends, generically, on the problem 
to be solved. Analyzing some propositions (3D harmonic oscillator and Coulomb 
wave functions) we employed the functions proposed by Fulcher |16| : these permit a 
simpler numerical treatment with less error possibilities when a numerical integration 
is eventually truncated 

That function described in expression (jl7p was also employed by other authors 
[TSl [T6l [TTI fT9] and it presents a similar analytical behavior to the Schrodinger wave 
function of the Coulomb potential: 



(17) 



where Lfc^''''^(/3r) are the associated Laguerre polynomials and the principal quantum 
number n of the bound state ri^^^^lj is written as n = k + 1. 

We are interested in obtaining a potential that describes, in a unique way, all the 
results of the interacting dynamics for charm systems as well as bottom systems, both 
with the same set of parameters. 

We introduce two more input parameters related to the electronic transition rates 
through the equation (fT8|) ,7]: 

r,-,+ (vi/) 4m,^{(3j ■ ^ 

The difference between charmonium and bottomonium resonances will depend on 
the parameters /3c and (3b in the wave functions. Experimentally we have |22j : 

T^-^+i^)) = 5.26 ± 0.37 keV, 

re-e+(T) = 1.32 ±0.05 keV, (19) 

indicating: 

2 

pc^pj^y. (20) 

Resuming the treatment, the calculations developed here are parameterized for: 
the charm and bottom quark effective masses {rUc and rr^^, the wave function 
parameter /3h and the potential parameters a, a, b and w [jjl. They are obtained 
through the above quoted mean square minimization. 

§ It provides a relative error of magnitude less than 10"^" when truncating the numerical integration. 
II The w of the exponential expression that parameterizes the interaction vector term may impose a 
Coulomb behavior to this term of potential, depending on the energy levels that we are computing. 
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Table 1. Phenomenological parameters for the charmonium 
unified effective model. 



rric 


= 1.108 GeV 


as 


= 0.739 


mt, 


= 4.727 GeV 


a 


= 0.0475 GeV^ 


13c 


= 0.148 GeV 


b 


= 0.088 GeV 


Pb 


= 0.390 GeV 


w 


= 0.550 GeV 



Table 2. Charmonium - ^ mass spectrum calculated in the use of a Coulomb 
plus Linear Potential. 



Effective Mass (MeV) 



niGson 




iOxp cr i mc ntr 


Theory 


D iffcrcncc 


*i 


1^50 


3096.9 


3096 


-1 




23^0 


3686.0 


3476 


-210 




3351 


3769.9 


3851 


+81 




4351 


4040 


4223 


+183 




5(?)35i 


4159 


4593 


+434 




6(?)35i 


4415 


4960 


+545 


Vc 


1^51 


2979.8 


3093 


+113 


XO 


l3Po 


3417.3 


3468 


+51 


XI 


l3Pi 


3510.5 


3468 


-42 


X2 


13P2 


3556.2 


3467 


-89 


7b 


l^Pl 


7 


3467 




7 


2^ Pi 


7 


3815 




? 


23 Po 


7 


3814 




? 


23 Pi 


7 


3815 




? 


23 P2 


7 


3815 




? 


33 Po 


7 


4160 




? 


33 Pi 


7 


4162 




? 


33 P2 


7 


4163 




? 


1^D2 


7 


3806 




? 


2^D2 


7 


4143 




? 


l^Di 


7 


3808 




? 


1^D2 


7 


3807 




? 


13 D3 


7 


3805 




? 


23Di 


7 


4145 




? 


23 D2 


7 


4145 




? 


23 D3 


7 


4143 





Measurement errors with last significant digit order |22) . 
? - Experimental data not accurately available. 



5. Results for mass spectrum of charmonia and bottomonia 

In accordance with the systeniatics described before, using the wave function proposed 
[16], we calculate the energy of rjc, ^'i, ^'2 and charmonium states and of Ti, T2, 
T3 and T4 bottomonium states, all with quantum numbers and experimental values 
(input data) described in tables and [31 We also use the values of expression as 
input data. 

With the expressions found we calculate the phenomenological parameters 
described in table [1] minimizing the function since we are interested in the 
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- bottomonium 
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Table 3. Bottomonium - T mass spectrum calculated in the use of a Coulomb 
plus Linear Potential. 



Effective Mass (MeV) 



meson 




Experiment 


Theory 


Difference 


Ti 




9460.4 


9608 


+149 


T2 




10023.3 


9931 


-92 


T3 


3^51 


10353.3 


10237 


— 116 


T4 




10580 


10533 


-47 


T5 




10865 


10825 


—40 


Tf"(?) 


6(?)3Si 


11019 


11113 


+96 


Xu) 


13P0 


9859.8 


9811 


-49 


Xbi 


T r-i 


9891.9 


9812 


-80 


Xb2 


T T-l 


9913.2 


9812 


-101 


^bO 


23 Po 


10231.1 


10042 


-189 


^bl 


23Pi 


10255.2 


10043 


—212 


/ 

Xb2 


23P2 


10268.5 


10044 


—224 


? 


l^So 


7 


9607 




? 


llPl 


7 


9812 




? 


2lPi 


7 


10043 




? 




7 


10270 




? 


33 Pi 


7 


10271 




? 


33 P2 


7 


10272 




? 




7 


9980 




? 


2^D2 


7 


10174 




? 


l^Di 


7 


9980 




? 


13D2 


7 


9980 




? 


13 D3 


7 


9980 




? 


23Di 


7 


10174 




? 


23 D2 


7 


10174 




7 


23 D3 


7 


10175 





parameter x' = 




(n is the number of input parameters). 



8 
i=l 

The term AM^ is the experimental error in the measurement of the mass Mi^^"^ . 
If AMi is less than 3 MeV , AMi — 3 MeV is adopted minimizing the influence of any 
input data [TT] , 

We found the minimum value x' — 42.6. 

In this way we obtain the parameters that characterize the quark - antiquark 
interactions in the charmonium - ^' and bottomonium - T resonances, obeying a 
Coulomb plus Linear behavior (see equation (fT5|) ). 

With those parameters we calculate the charmonium and bottomonium mass 
spectrum for several quantum numbers. Some authors [H O [17] suggest that the 
calculation of transition rates should be done in parallel to mass spectroscopy but the 
results, confirmed by ours, corroborate the fact that the mass spectroscopy is not the 
best way to do so. The mass spectroscopy results are described in tables [2] and [3l 

With these values, considering just the data with quantum numbers exactly 
confirmed through the experimental measurements, we observe a mean deviation of 
3.54% in the charmonium mass spectrum and a mean deviation of 1.14% in the 
bottomonium mass spectrum, when we compare with the experimental data. 



the 



AM, 



(21) 
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Since we were working all the time with a reduced and unified number of 
parameters for charm and bottom families, the results obtained describe the mass 
spectroscopy in a way as good as other results [SI [HI [HI [Sj obtained with a more 
complex and extensive analysis. 

Other aspects also deserve some emphasis. It is clear, with corroboration of some 
results already obtained, that the spectroscopy study with different sets of parameters 
for each meson family would provide better results despite the fact that experimental 
data with rigorously accurate quantum numbers are not sufficient to establish a trusty 
statistics. 

Attending to the accuracy of the nonrelativistic approximation adopted for 
the kinetic term, it is quite natural that we observe a minor mean deviation for 
bottomonium states since charmonium states could be considered semirelativistic. 

Another important observation is that, like results of other authors [II [6l [8l [HI 
[SJ [S] , there is an experimental disagreement relative to degenerancy in total angular 
momentum J. Observing the results in tables [D and [31 we can conclude that Breit - 
Fermi interaction terms are not sufficiently accurate to describe states with different 
values of J (J — I — 1, I and I + 1). 

6. Conclusion 

When a study of mass spectroscopy of quark bound states is developed, the first 
aspect that will establish its success is obviously the correspondence of the theoretical 
model to the experimental measurements. In this way, the potential model developed 
and applied here permits a broad and accurate description of the mass spectrum of 
charmonium and bottomonium systems when we attempt to the mean deviation of all 
resonances. 

Adopting an effective Hamiltonian formalism built on a restrict set of eight 
parameters!^ proved to be a significant simplification in the calculations here 
developed. Other authors [H [3 [9l [11] employ more complex developments that 
often use more than twenty parameters. In spite of the use of some complex 
developments, none of them is capable to describe states with different values of J 
(J — I and I + 1) accurately. Some of them [8l [151 [12] do not even bother about 

such calculations involving the total angular momentum quantum number. Others 
[51 [S] [TT] obtain results with deviations as large as those presented in this paper. 

The first simplification occurs because some nonlocal terms that would appear 
in the explicit potential could be "absorved" by the linear parameters. The mass 
spectroscopy of hadrons following a potential model is not the best way to indicate 
the exact form of the nonlocal terms. 

Searching for a simpler model, we succeed in describing the kinetic term in a 
nonrelativistic approximation in which, differently to other models [6l 18] fT5 ] fT6 l I17j. 
we adopt a kinetic term accurate in the energy region below 10 GeV. It minimizes the 
difference to the exact kinetic term. 

The advantage of such procedure is that the obtained values for po^ are 
independent of any energy state used as input data. In this way the extension of 
using the same value of po^ for any other resonances of that particular bound state 
is immediate. In comparison, Jaczko and Durand do not extend their calculation in 

1 We remind that fSc and (3), are linked by expression l|20|l . so there are only seven independent 
parameters. 
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|17| beyond the four fower energy states which otherwise were extensively studied in 
their paper. The quoted authors [17] had chosen to keep the kinetic term depending 
on one more parameter (e') in the effective Hamiltonian. 

AU that simphfication together with a reduced number of eight parameters aUows 
for the use of an also reduced number of input data. With this, the higher energy 
levels do not present a direct constraint with the experimental data, a fact that would 
force a better agreement. 

It should be emphasized that the set of parameters is common to both families 
of charmonia and bottomonia. Lower deviations from the experimental values could 
have been achieved if a different set of parameters had been adopted for each family 
separately. 

Despite the relative success, some questions concerning potential models yet 
remain unanswered: Why do we not try to solve directly (numerically) a relativistic 
bound state equation instead of following an effective model? 

Generically, we can conclude that there is a great number of possibilities to 
solve a quark bound state problem, sometimes through analytically or numerically 
complicated calculations, but none of them is exact. We hope that some direct 
or indirect results of QCD can provide new guides to produce more accurate 
phenomenological models. 
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